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Input file Fbh62112FL.seq; Output File 62112. trans 
Sequence length 2452 

M SO 3 

aCTCTCTCTG ^^ ATC AGC GGC 9 

CGLFLRTTAAARACRGLVVS 23 
TGC QGG CTC TTC CTG CGC ACC AOS GCT GCG GCT CGT GCC TGC CGG GGT CTG GTG CTC TCT 69 

TANRRLLRTSPPVRAFAKEL 43 
AOC GCG AAC OQG OQG CEA CTG OQC ADC AGC OOG OCT CTA CGA OCT TTC GCC AAA GAG CTT 129 

FLGKIKKKEVFPFPEVSQDE 63 
TTC CTA GGC AAA ATC AAG AAG AAA GAA GTT TTC CCA TTT CCA GAA GTT AGC CAA GAT GAA 189 

LNEINQFLGPVEKFFTEEVD 83 
CTT AAT GAA ATC AAT CAG TTC TTG GGA CCC GTG GAA AAA TTC TTC ACT GAA GAG GTG GAC 249 

SRKIDQEGKIPDETLEKLKS 103 
TCC OGA AAA ATT GAC CAG GAA GOG AAA ATC CCA GAT GAA ACT TTG GAG AAA TTG AAG AGC 309 

LGLFGLQVPEEY GGLGFSNT 123 
CTA OQG CTT TTT GOG CTG CAA GTC CCA GAA GAA TAT GGT GGC CTG GGC TTC TCC AAC ACC 369 

MYSRLGEIISMDGSITVTLA 143 
ATG TAC TCA AGA CTA GGG GAG ATC ATC AGC ATG GAT GOG TCC ATC ACT GTG AOC CTG GCA 429 

AHQAIGLKGIILAGTEEQKA 163 
GCG CAC CAG GCT ATT GGC CTC AAG QGG ATC ATC TTG GCT GGC ACT GAG GAG CAG AAA GCC 489 

KYLP KLASGEHIAAFCLTEP183 
AAA TAC TTG OCT AAA CTG GCG TCC GGG GAG CAC ATT GCA GCC TTC TGC CTC ACG GAG OCA 549 

ASGSDAASIRSRATLSEDKK 203 
QCC AGT GGG AGC GAT GCA GCC TCA ATC CGG AGC AGA GCC ACA CTA ACT GAA GAC AAG AAG 609 

HYILNGSKVWITNGGLANI F 223 
CAC TAC ATC CTC AAT GGC TCC AAG GTC TGG ATT ACT AAT GGA GGA CTG GCC AAT ATT TTT 669 

TVFAKTEVVDSDGSVKDKIT 243 
ACT GTG TTT GCA AAG ACT GAG GTC GTT GAT TCT GAT GGA TCA GTG AAA GAC AAA ATC ACA 729 

AFIVERDFGGVTNGKPEDKL 263 
GCA TTC ATA GTA GAA AGA GAC TTT GGT GGA CTC ACT AAT GGG AAA CCC GAA GAT AAA TTA 789 

GIRGSNTCEVHFENTKIPVE 283 
GGC ATT CGG GGC TCC AAC ACT TCT GAA GTC CAT TTT GAA AAC AOC AAG ATA OCT GTG GAA 849 

NI LGEVGDG FKVAMNI LNS G 303 
AAC ATC CTT GGA GAG GTC GGA GAT GGG TTT AAG GTG GCC ATG AAC ATC CTC AAC AGC GGC 909 

RFSMGSVVAGLLKRLIEMTA 323 
CGG TTC AGC ATG GGC AGC GTC GTG GCT GGG CTG CTC AAG AGA TTG ATT GAA ATG ACT GCT 969 

EYACTRKQFNKRLSEFGLIQ 343 
GAG TAC GCC TGC ACA AGG AAA CAG TTT AAC AAG AGG CTC AGT GAA TTT GGA TTG ATT CAG 1029 

EKFALMAQKAYVME SMTYLT 363 
GAG AAA TTT GCA CTG ATG GCT CAG AAG GCT TAC CTC ATG GAG AGT ATG ACC TAC CTC ACA 1089 

AGMLDQPGFPDCSIEA AMVK 383 
GCA GGG ATG CTG GAC CAA CCT GGC TTT CCC GAC TGC TCC ATC GAG GCA GCC ATG GTG AAG 1149 

Fig. 1A 
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Protein Family / Domain Matches, HMMer version 2 



Searching for complete domains in PFAM 

hmrpfam - search a single seq against HMM database 

HMMER 2.1.1 (Dec 1998) 

Copyright (C) 1992-1996 Washington University School of Medicine 
HMMER is freely distributed under the GNU General Public License (GPL) 



HMM file: 
Sequence file: 



/prod/ddm/seqanal/PFAM/pf am6 . 4/Pf am 
/procVddm/wspace/orf anal/oa- script . 26629 . seq 



Query: 62112 

Scores for sequence family classification (score includes all dcmains) : 
Model Description Score E -value N 



Acyl -CoA_dh_M 
Acyl-CoAdh 
Acyl - CoA dh N 
Polysacdeacet 



Acyl-OoA dehydrogenase, middle dcmain 
Acyl-CoA dehydrogenase, C- terminal dcm 
Acyl-OoA dehydrogenase, N- terminal dcm 
Polysaccharide deacetylase 



153.0 
152.1 
73.7 
-43.7 



3.5e-42 
9.6e-42 
4.2e-19 
1.8 



1 
1 
1 
1 



Parsed for dcmains 
Model 



main 


seq-f 


seq-t ' 


hmm-f 


hmm- 1 




score 


E-value 


1/1 


85 


177 


29 


132 


.] 


73.7 


4.2e-19 


1/1 


179 


2o6 • • 


1 


106 


[] 


153.0 


3.5e-42 


1/1 


290 


441 . . 


1 


156 


[] 


152.1 


9.6e-42 


1/1 


432 


580 . . 


1 


150 


[] 


-43.7 


1.8 



Acyl - CoA_dh_N 
Acyl -CcAdhM 
Acyl-CoAdh 
Polysacdeace t 

Alignments of top- scoring dcmains: 

Acyl -CoA dh_N: dcmain 1 of 1, frcm 85 to 177: score 73.7, E = 4.2e-19 

* - >RRvDksgef PalrelikaLgqlGllginvPEeyGGaGad . . ylaRRn 
R++D++g+ P e +++L lGl4g+ vPEeyGG+G +++ ++ 
62112 85 PKIDQEGKIP- - DETLEKLKSDGLFGIjQW 126 



62112 127 



LHAQVaalviEElarvcAstgv^ 

+ E+ ++s v+l++h ++g+ +i+ +G+eEQK+kyIp+ 
RliG^ISMDGSriVIljAAHQ^ 168 



62112 



ltsGdligao* 
l+sG++i+a 
169 LASGEMAA 



177 



Acyl-CoA dh_M: dcmain 1 of 1, frcm 179 to 286: score 153.0, E = 3.5e-42 

*->AlTEPgAGSDvgSlkTtAekkEGd. . dyil^lSsKmWITOGgqAdwyi 
+1TEP +GSD++S++ +A+ d+++yiI25GsK+WITOQg A++++ 
62112 179 CLTEPASGSDAAS IRSRATLS - EDkkHYTI2*3SKVWITWQGLANI FT 224 



62112 



VlAvT. . . Dpakkvpgkkgi t^lVekdt^f siGkKedKDGlRgS<frcE 
V+A+T+ D + + k itaF+Ve+d+ G++ Gk+edKLG+RgS+TcE 
225 VFAKTewDSDG- - SVKDKTTAF1 Vh^FQGVTtK^EDKLGlRGSZTTCE 272 



62112 



LiFEDvrvPesniL< - * 
+ FE+ ++P +niL 
273 VHFENIKIPVENIL 



286 



Fig. 2A 
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Acyl-OoA_dh: dcmain 1 of 1, fran 290 to 441: score 152.1, E = 9.6e-42 

G+<3Fk+anH- L+ +R+ +++ G1+ + ++ +++Ya RkqF k+1 
62112 290 <3*3FKVAMOTI2JSGRF 336 

adfQliQfkLAdMatkLE^aRllvYr . GedAKEAIptskeaam 
+f liQ+k+A Ma k ++ +++Y +a d+-K3 + ++s eaam 
62112 337 SEFXaLIQEKEMWAQKAYV DCSIEAAM 381 

AKlfaseaAinqvatudAvQilGGvGYtJcd^ 

+K f +seaA q +++A+QilQG GYt dyP eR +RD +i I BCT+EI 
62112 382 VKVFSSEAAW2C^SEAI^ 431 

qrlvIaRall<-* 
r la + 1 
62112 432 LRMYIALTGL 441 

Polysacdeacet: ckxna±n 1 of 1, frcm 432 to 580: score -43.7, E = 1.8 

* - >ddksvyLTFDD^^ 

+++++LT + ++ + T+r+ + Lk+ kv + G++ +d 
62112 432 LRMYIALTGU3HAG-- -RILTTRI-HEL^^ - 473 

dlarrivkeGHei^SHtwsHPdlt tl 

+ r v+ G gNH+ HP 1+++ ++ ++++ +++ ++ + + 
62112 474 -SI/3RTVDIX5-LTG^GWH 521 

toeqirdeiertneaiiqatggatptlfl^^ 

t +++ + r +++++i+++g t++l R+ s+s ++G1+ 

62112 522 TmE^IaVIXKVim SRSIRIGLRN 560 

vlWdvDprEWsvragadaivdavlqaa< - * 
+ D v ++ v a+lq+ 

62112 561 H EHEVIiLANTFCVEAYLONL 580 



Fig. 2B 
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Protein Family / Domain Matches, HMMer version 2 
Searching for complete domains in PFAM 
hnnpfam - search a single seq against HMM database 
HMMER 2.1.1 (Dec 1998) 

Copyright (C) 1992-1998 Washington University School of Medicine 
HMMER is freely distributed under the GNU General Public License (GPL) . 

HMM file: /prod/ddm/seqanal/PFi^pfam5.4/Pfam 
Sequence file : /prod/ddm/wspace/orf anal/oa-script . 17193 . seq 

Query: 62112 

Scores for sequence family classification (score includes all domains) : 
^^l. _ - Description Score E-value 



Acyl-CoA dh Acyl-CoA dehydrogenase 399.8 1.8e-116 

Polysac_deacet Polysaccharide deacetylase -43.7 l 

Parsed for domains: 

Domain seq-f seq-t hmm-f hmn-t score E-value 

Acyl-CoA_dh 1/1 85 438 .. 29 394 .] 399.8 1.8e-116 

Polysac deacet 1/1 432 580 . . 1 150 [] -43.7 l 

Alignments of top- scoring domains: 

Acyl-CoA dh: domain 1 of 1, from 85 to 438: score 399.8, E = 1.8e-116 

* - >RRvDksgefPlrelikaLgklGllginvPEeyGGaGad . . ylaRFmL 
R++D++g+ P e +++L lGl+g+ vPEeyGG+G +++ ++ 
62112 85 RKIDQEGKIP -DETLEKIJCSLGLPGLQ\^EErYQGLGFSntMyS 126 

HAQVaal viEEl arvcAs tgvl 1 svhs sLgqnpi 1 r f GseBQkkkyLpql 
+ E+ ++s v+l++h ++g+ i+ +G+eEQk+kyLp+l 
62112 127 RI/3EI ISMDGS ITVTLAAHQAIGLKGI ILAG7TEEQKAKYLPKL 169 



Fig. 2C 
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t sGdl igafAlTEPgAGSDvgS ikTtAekkEGd . . dyaI2^sKirtVIlNGg 
+sG++i+af+lTEP +GSD++Si+ +A+ d+++yiI^GsK+WITOGg 
62112 170 ASGEHIAAPCLTEPASGSDAASIRSRATLS - EDkkHYII2$GSKVWITOB9G 218 

qAdwyiVlAvT. . . DpakkvpgkkgitaFlVekdtpGf siGKKedKDGlR 
A++++V+A+T+ D + + k itaF+Ve+d+ G++ Gk+edKLG+R 
62112 219 LANIFTVFAKTevvDSDG- - SVKDKITAFTVEE^PGGVTOGKPEDKI/SIR 266 

gS<±TcELiFEE>vTrvPesniljGeBGeGFk^ 

gS_TcE+ FE+ ++P +niLGe G+GFk+aM+ L+ +R+ +++ G++ 
62112 267 GSNTC£WHFEOTKIPVEtn:i/3^^ 316 

gAldeAiiiYAkqRkqFGkplae f Ql iQ fkLAdMAtkLEaaRl IvYr aAwl 
+ ++ +++YA RkqF k+1 ef liQ+k+A MA k ++ +++Y +A 
62112 317 RLIEMIAEYACTRKQFNKRIi^ 366 

adr . GedAKEAI^tskeAAMAKlfAseiArnkvatdAvQil(^^ 
d++G + ++s eAAM+K f +se+A + +++A+QilGG GYt dyP 
62112 367 LDQpGFP rX^IEAAMVKVFSSEAAWQCVSE^ 411 

veRfyRDAki tqlYECTsEIQr 1 vIaR< - * 
eR +RD +i I BGT+EI r la 
62112 412 YKPTT .ftytpttx .tttpyttntktt .pmytal 438 

Polysac deacet: domain 1 of 1, from 432 to 580: score -43.7, E = 1 

* - >ddksvyLTFDDGPnAApayTprl 

+++++LT + ++ + T+r+ + Lk+ kv + G++ +d 
62112 432 LRMYTALTCLQHAG- -RILTTRI -HEIJCQAKVSTVMDTVGRRIjRD- - 473 

dlarrivkeGHei^IHtwsHPdlt tl 

+ r v+ G gNH+ HP 1+++ ++ ++++ +++ ++ + + 
62112 474 - SDGRTVDLG - LTGNHGWHPSLAdsankf een tyc f gr tve 1 1 1 1 r f GK 521 

teeqirdeiertaeaiiqatggatptlfRpPYGewsetvlsasaklGlta 
t +++ + r +++++i+++g t++l R+ s+s ++G1+ 

62112 522 TIMEBQLVIiKRVANILINLYG - MTAVLSRA SRSIRIGLRN 560 

vlWdvDprDWsvragadaivdavlqaa< - * 
+ D v ++ v a+lq+ 

62112 561 H DHEVIiANTPC^/EAYljQNL 580 



Fig. 2D 
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